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The acoustic emission from breaking a bamboo chopstick or a bundle of spaghetti is found to 
exhibit similar behavior as the famous seismic laws of Gutenberg-Richter, Omori, and Bath. By use 
of a force-sensing detector, we establish a positive correlation between the statistics of sound intensity 
and the magnitude of tremor. We also manage to derive these laws analytically without invoking the 
concept of phase transition, self-organized criticality, or fractal. Our model is deterministic and relies 
on the existence of a structured cross section, either fibrous or layered. This success at explaining 
the power-law behavior supports the proposal that geometry is sometimes more important than 
mechanics. 

PACS numbers: 62.20.mt, 91.30.Px, 89.20.-a, 43.50.-|-y 


Fracture [T] is a complex phenomenon with many 
interesting properties, such as crack propagation^], 
breakdown [3|, and self-affine fractals in the crack surface 
morphology [3]. This kind of studies started in seismol¬ 
ogy and attracted the attention of physicists and mate¬ 
rial scientists. Crackling noiseO Ej in the acoustic emis¬ 
sion from fracturing wood plate [7], paper [8l[9|, rock[T0|, 
concrete[llj. and charcoal [T^ has been found to exhibit 
similar statistical properties as earthquake [TOl IHj. On 
the other hand, the slow growth of a single crack in a 
fibrous sheet has been studied by the lattice model[9] 
which incorporates thermodynamics to describe the tem¬ 
perature dependence. To understand the failure process, 
the fiber bundle model[T5| has also been extendedpffi] 
by introducing time-dependent damage accumulation of 
fibers to capture the stochastic nature of fracture. 

In contrast to 2- and 3-D samples in previous frac¬ 
ture experiments, we shall adopt the 1-D like chop¬ 
stick because it rarely fractures in multiple places at the 
same time, which makes the interpretation of data eas¬ 
ier. Chopsticks are easily accessible because they are 
indispensable implements for most Asian families and 
carry special cultural meanings |17j in Chinese tradition. 
They are usually made by bamboo in Taiwan and China. 
Bamboo exhibits a fibrous cross section. Its excellent 
compressive strength, endowed from the vascular bun¬ 
dles scattered throughout the stem in the cross section 
instead of in a cylindrical arrangement, enables us to col¬ 
lect about 400 crackling sounds from breaking one chop¬ 
stick, as opposed to just one sound from a chalk or twig. 

Our main sample is the cheap and omnipresent bam¬ 
boo chopstick of diameter 0.67cm and length 20.5cm. 
Sticks are also taken directly from the bamboo tree to 
make sure our conclusions are independent of the length, 
freshness, and processing to turn them into chopsticks. 
A bundle of spaghetti m is also used to check the gen¬ 
erality of our model because it shares similar cross sec¬ 
tion as bamboo (see Fig0. The fracture machine in 
Figj^ is placed inside a soundproof chamber with foam 
rubber plank on the interior to avoid echo. A Sony 
ECM166BC microphone that is connected to a Sony 
ICD-PX333 recorder picks up fracture sound. By tight¬ 



FIG. 1: (color online) Panels (a) and (c) show the similar 
structure in cross section between a bundle of spaghetti and 
a bamboo chopstick. Panels (b) and (d) are the side view 
during fracture. 


ening the screw, the movable metal block was driven and 
compressed the chopstick. The bolted plate is to secure 
the bent chopstick from slipping and flying off. The chop¬ 
stick was flipped into a horizontal position as soon as it 
became bent to prevent fracturing near the plate edge. 
Crackling sound was recorded at a sample rate of 44100 
points per second in 16-bit precision. The amplitude was 
measured in computer units and the maximum amplitude 
Amax = 2^^ — 1. The gain of sound card was constant and 
samples were positioned at a distance of 5 cm from the 
microphone. We kept the process of fracture at a steady 
velocity of about 200 s per chopstick. 

Average amplitude of background noise was 3 x 10“^ 
as normalized by Amax) and thrice this amount was set as 
the noise threshold. The c code algorithm automatically 
integrated the sound intensity every 200/44100 s. When 
the value exceeded background noise, beginning of a new 
pulse was marked, as in Figj^c). Whenever a dilemma 
arose at distinguishing a long pulse from two overlapping 
pulses as Figj^d), we resorted to a smaller time step to 
examine the intermediate area by including just six am¬ 
plitude peaks of the oscillation within the pulse(s). Since 
this value is expected to decrease as a pulse fades, a sud¬ 
den switch to increasing function indicates the beginning 
of a second pulse. 
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FIG. 2: (color online) Panels (a,b) show the fracture machine. 
FSR is the acronym for force-sensing resistor, and GND is the 
shorthand for ground. The blue dotted arrow represents the 
arm of force, and the red solid arrow the force exerted on one 
end of the chopstick. Panel (c) shows a typical sound pulse, 
and (d) an example when one long pulse or two overlapping 
pulses need to be distinguished. 


After squaring the amplitude and integrating over the 
duration time to obtain the energy for each sound pulse, 
we grouped the data into a histogram of probability func¬ 
tion versus pulse energy. A full-logarithm was then taken 
in Figj^a) and (d) to reveal a good alignment, indicative 
of a power law: P{E) ^ IjE^ that is reminiscent of the 
Gutenberg-Richter law[Tn]. The exponent, b, for chop- 
stick is about 1.45 by use of the maximum likelihood 
estimation (MLE). When we sliced the chopstick hori¬ 
zontally into two halves, the power law persisted but the 
upper (lower) hemisphere became b = 1.48 (1.40). The 
empirical value of b for earthquake is commonly close 
to 1.0 in seismically active regions, but may vary in the 
range of 0.5 and 2 depending on the source environment 

of the region [20] • 

The aftershock rate in Fig|^b) and (e) was defined 
as the number of aftershocks divided by the time win¬ 
dow of 100ms, against the elapsed time, tEiapsed- And 
aftershocks referred to the sound pulses on the time se¬ 
ries between two local maxima in pulse energy; i.e., the 
main shocks. Since the number of local maxima roughly 
matches that of vascular bundles, we believe the fracture 
of each bundle gives rise to one main shock and a few 
aftershocks. A shifted power law, rate= fc/(c-I-tsiapsed)^ 
where k and c are constants, similar to the Omori’s 
law|2I] was deduced. The exponent in Figj^b) was de¬ 
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FIG. 3: (color online) Panels (a, b, c) are for crackling sound 
from the bamboo chopstick, while (d, e, f) from a bundle of 
spaghetti. They mimic the Gutenberg-Richter law, Omori’s 
law, and Bath’s law for earthquake. The power-law exponent 
in (a) and (d) is determined to be 1.45 and 1.35. The expo¬ 
nent and shift in (b) and (e) are {p = 1.68, c = 0.07) and 
(3.53, 0.027). The relative magnitude M in panels (c) and (f) 
is defined as the average of log(i?ms) — log(ii'ia) where Ems 
and Eia, denote respectively the energy of main shock and its 
largest aftershock. Error bars represent standard deviations. 


termined by MLE to be p = 1.68, slightly larger than 
the region of 0.75-1.5 set for earthquakes. Finally, after 
an initial dip in Figj^c) and (f) that is known to ex¬ 
ist in real datajlH I23] 7 the M-value became independent 
of the main-shock magnitude and saturated at 1.7-1.8 in 
Fi j^c) - somewhat higher than 1.1-1.2 for the Bath’s 

In spite of a multitude of workers who have compared 
the acoustic emission from laboratory material to the 
shock energy of earthquakes, the reasonable assumption 
that these two energies are correlated was rarely proven. 
We thus decided to afRx a force-sensing resistor (FSR) 
to the end of chopstick in Figj^and compared the timing 
and statistical behavior of the discrete events of tremor 
and acoustic emission. When the reaction force by the 
chopstick increased, the resistance of FSR (Interlink Elec¬ 
tronics FSR402) decreased from an initial value higher 
than lOM ohm and generated a nonzero voltage, 14- An 
Arduino UNO board then collected the measurement of 
14 at a sample rate of 100 points per second. 

Any fracture is expected to cause a sudden drop in re¬ 
action force against the FSR, but not necessarily accom¬ 
panied by a detectable sound. Therefore, there is no one- 
to-one correspondence between these two events, with 
there being more tremors than sounds. Any discussion 
of their correlation has to be effectuated through compar¬ 
ing their statistics. We thus grouped the nonzero values 
of the time derivative of force, dF/dt, into a histogram. A 
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power-law distribution P(dF/dt) ^ l/{dF/dt)^ emerged 
with /? = 1.31 in Fig|^ This is to be contrasted to 
Figj^a). The fact that both figures exhibit power-law 
behavior with similar exponent and valid range of pa¬ 
rameter confirms a positive correlation between energies 
of the quake and the crackling sound. The chaotic dis¬ 
tribution when dF/dt exceeds lON/s was caused by the 
drop of sensitivity in FSR when the force change is larger 
than ION. Limited by the sensitivity of FSR, we believe 
the number of small events it picked up was underesti¬ 
mated which probably explains why /3 = 1.31 is slightly 
smaller than b = 1.45 in Figj^a). 



dF/dt (N/s) 

FIG. 4: (color online) Full-log plot of the number of events 
versus the magnitude of force change picked up by FSR. Data 
were fit to a power law (the red dotted line) with an exponent 
roughly equal to that of crackling sound in Fig|^a). 

A minimal model will be presented to explain qual¬ 
itatively why the major seismic laws should appear in 
the humble bamboo chopstick. In contrast to numeri¬ 
cal simulations that are popular for the study of fracture 
in disordered media[2S], we believe analytical derivations 
can provide more insights. Since the Gutenberg-Richter 
and Bath laws concern only the energetics, they will be 
tackled separately from the Omori’s which requires dy¬ 
namics. For simplicity, we only consider a rectangular 
cross section, but other shapes can be easily generalized. 

Two features in Figj^merit attention. First is the pro¬ 
gressive shortening of effective length, - the length of 
beam neutral axis that is actually bent. Second, the end 
segments of the chopstick remain straight because they 
are thicker and less susceptible to bending. To nail down 
the relevant factors and enable analytic solutions, we will 
focus on the bending energy in the static beam equation, 
which includes tensional and compressional stresses in 
the cross section external and internal to the beam neu¬ 
tral axis, respectively: 

(s/) 

where ATb ~ Y{na)^ and Y denotes the Young’s modu¬ 
lus. As exemplified by Figj^b, c), the fractured segment 
of bamboo remains rather straight, which demonstrates 


that the plastic region preceding fracture is very narrow 
and the elastic form of Eq.Q can be extended to large 
strains. As fracture proceeds in the tensional region, n 
decreases and In is shortened progressively as evidenced 
by Figga, b, c) and our YouTube clip|26j. 



FIG. 5: (color online) Yellow double-arrowed lines in panels 
(a,b,c) highlight the shortening of effective length. In panel 

(d) , fiber on the n-th layer is stretched and on the verge of 
fracturing. The Rn represents its radius of curvature, and 6n 
the angle spanned by the bent section (in dotted red line). 
Note the end segments (in blue dashed line) remain straight, 
resembling the V-shape configuration in Ref. m- In panel 

(e) , the reduction ratio of effective length is determined for 
different time stages. 


Geometric relations can be written down from Figj^d) 


F — (-^n F Ti(Fj On 
In = [Ru F y) On 


Subtracting them gives 


^ n _ a; - w 

2 ’ 


( 2 ) 


(3) 


where the threshold strain for fracture, e = Ain/In, is a 
material property. Simple algebra gives 


In. 


na 
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since Rn S> na. For a rectangular cross section, the y- 
integration simply gives a constant and the energy stored 
in the n-th layer is roughly 


En 
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(5) 
































4 


Since the stress is larger away from the beam neutral 
axis, the equipartition approximation of En/Etot ~ 1 /n 
is an underestimation. 

The fracture process was divided into four time stages 
in Figj^e). If there are 400 pulses, take the ratio of 
hoi/hooj where /301 and ^400 represent the effective length 
measured upon the 100 th and the first pulses, and dehne 
it as where A is the mean value. Figurej^e) shows 
that A is roughly constant. The deviation of A value and 
a large error bar in the final stage is likely caused by the 
uncertainty in pulse number because crackling sound be¬ 
comes very feeble. From these observations, the effective 
length can be assumed to obey 

ln=l- ( 6 ) 

where I denotes the initial length and A < 1 is a material 
property. Similar concept has been used to explain the 
power-law behavior of the crackling sound from a crum¬ 
pled thin sheet [28]. 

The total number of crackling sound equals 

Divide both sides by S^t unity on the left side 

to conserve probability while the occurrence rate can be 
identified as P{En) oc dn/dEn ^ l/En by use of Eqs.([^ 
1^. This reproduces the power law in Figj^a). If we 
increase the ratio En/Etot = ^/n, the exponent will be¬ 
come larger and closer to the empirical value b = 1.45 in 
Figj^a). When the rectangular cross section is changed 
to, say, southern hemisphere of a horizontally halved 
chopstick, the width and hber number of each layer de¬ 
crease as fracture propagates downward from the top. 
Since top layers emit stronger pulses according to Eqs.([^ 
1 ^, there is now more count of loud sound than weak 
sound. And we expect Eig.3(a) to level off and give a 
smaller b. This is consistent with our observation. 

To derive the Omori’s law, we need the time interval 
between successive aftershocks to calculate the elapsed 
time in Fig.3(b). By assuming a constant angular accel¬ 
eration during this short interval, we can estimate this 
time, Atn for the torque to render a displacement of A0„: 

A9n = On-i -0n = ^ (FZ sin 6 >„//)(At„)^ ( 8 ) 

where I denotes the moment of inertia. Since the bending 
force, F, comes from far ends of chopstick, the arm of 
force equals I, not In- The sin0„ factor is due to the fact 


that the force is not perpendicular to its arm. Finally, 
since is mostly small, sin0„ ~ On- As n ^ 1, we can 
replace A0„ by dOn/dn and obtain from Eq.([^: 

(At„)^ oc In A -I-- (9) 

n 

by use of Eqs.(|^ |^. This shows that the time inter¬ 
val to the next aftershock increases with the number of 
broken hbers, N — n- A simple picture to understand 
this slowdown in braking is that the tensional stress ex¬ 
ternal to the beam neutral axis is lessened, as the chop¬ 
stick gets thinner. Consequently, it takes a larger angu¬ 
lar displacement and longer time to reach the threshold 
strain for fracture. The aftershock rate can be calcu¬ 
lated by dividing the number of aftershocks, rims ~ 
by the time lapse from their corresponding main shock, 
Atn^,-i + Atn ^,-2 + ■ • • + Atn where rims > 1 labels 
any one of the main shocks. Since the maximum number 
of aftershocks, rims ~ is in the ballpark of 8 , we can 
Taylor-expand the denominator by the small parameter, 
1 — {n/rims) ^ 1 and obtain 


rate = 


’^ms ^ 


'^ms 

D' + {D"/2)[l-{n/nms)] 


( 10 ) 


where the D' and D" denotes the first and second deriva¬ 
tive. Since rims — n is roughly a measure of time, Eq.( 10) 
reproduces the shift and power law of Figj^b). 


Since F„/F„_i ~ 1/A is independent of the value of n 
and insensitive to the approximation of En/Etot = ^/n^ 
the Bath law is naturally guaranteed. 

In conclusion, we establish a complete parallel between 
the crackling noise of a common bamboo chopstick and 
a bundle of spaghetti and the fundamental seismic laws. 
The statistics of these acoustic events is shown to corre¬ 
late with that of tremor. The htting function and corre¬ 
sponding parameters are determined by the rigorous sta¬ 
tistical method of Akaike Information Criterion [28l [29] 
and MLE. We succeed at deriving these laws analytically 
without invoking the concept of phase transitionjS], self- 
organized criticality jSni I3I|, or fractaljUj. Our models 
are based mainly on a structured cross section, which 
can be either hbrous or layered. This shift of empha¬ 
sis from mechanics to geometry to explain the power-law 
behavior is in conformance with the proposal of Ref. [32] . 
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